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«) , Abstract Solar flare accelerated electrons escaping into the interplanetary space 

■ and seen as type III solar radio bursts are often detected near the Earth. Using 
^ . numerical simulations we consider the evolution of energetic electron spectrum 

■ in the inner heliosphere and near the Earth. The role of Langmuir wave gener- 
, ation, heliosphcric plasma density fluctuations, and expansion of magnetic field 

lines on the electron peak flux and fluence spectra is studied to predict the 
\ electron properties as could be observed by Solar Orbiter and Solar Probe Plus. 

Considering various energy loss mechanisms we show that the substantial part 
of the initial energetic electron energy is lost via wave-plasma processes due to 
plasma inhomogeneity. For the parameters adopted, the results show that the 
electron spectra changes mostly at the distances before ~ 20 Rq. Further into 
the heliosphere, the electron flux spectra of electrons forms a broken power-law 
relatively similar to what is observed at 1 AU. 

^ ■ Keywords: Energetic Particles, Propagation; Flares, Dynamics; Radio Bursts, 

' Type III; Solar Wind; Waves, Plasma 

1. Introduction 

C^ - . Solar flares can efficiently accelerate a large number of electrons to sub-relativistio| 

energies. These energetic electrons are prominently observed via their hard X-ray 
emission (see e.g. lHolman et 001201 11 iKontar et «UI2011l as recent reviews), but 
can also escape into interplanetary space. Energetic electrons that escape flaring 
regions can be detected in-situ by near-Earth particle detectors ( |Lin, 20l"Tj ). 
1 Detection of escaping energetic electrons closer to the Sun is via their char- 

acteristic emission of both solar and interplanetary type III radio bursts. As 
the electrons travel further away from the Sun the faster electrons overtake the 
slower ones creating a positive gradient in velocity space. When the number of 
1 energetic electrons is sufficiently high Langmuir waves can be generated. The 
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presence of a high level of Langmuir waves can be viewed via the associated 
plasma emission, seen as type III radio bursts. 

For the electrons exciting type III bursts the transport is governed chiefly by 
beam-plasma interactions involving electron-Langmuir wave interactions ( Zhclcznyakov and Zaitsev, 1970 | 



Lin et al, 1981 Melrose, 1990). This process is fast and the characteristic length 



of interaction or Langmuir wave generation by electrons with velocity v and 
electron density n(v) scales as ~ vn p /[u)p e n{v)\, where u> pe is the electron plasma 
frequency and n p is the plasma density. The complicating aspect of the electron 
transport and Langmuir wave interaction is that Langmuir waves are effectively 
scattered and refracted by plasma density fluctuations. This results in a fast 
change of Langmuir wave spectrum, which in turn affects the overall evolution 
of the electron stream travelling from the Sun to the Earth. Historically, sim- 
ulations have been used to model the processes of electron beam generation of 
Langmuir waves (e.g. ITakakura and Shi bahashi, I1976P as well the role of den- 
sity inhomogeneities (e.g. Nishika wa and Ryutov| 119761 IGoldman and Dubois! 
115551 IKontarl I2001bl |Li, Robinson, and Cairns} I2D051 IReid and KolrtaTl I50TU1 
IZiebell et aU 12011^ . 

Recently IReid and Kontarl (|2010p . IKontar and Reidl (|2009p have shown that 
because of the inhomogeneous plasma of the solar wind, the spectral break in the 
peak flux and fluence appears as the electrons travel from the Sun to the Earth. 
As a result, the electron flux spectrum (peak and fluence) at 1 AU appears close 
to a broken power-law with a typical spectral index below the break around —2. 
Importantly, the exact value is dependent on the level of density fluctuations, so 
that the spectral index of electrons below the break is higher for a higher level of 
density fluctuations. Above the break the electron spectrum is weakly affected by 
the Langmuir waves, so can be considered as 'scatter-free' transport. However, 
these electrons are scatter-free only in respect to Langmuir waves and could be 
affected by other plasma waves (e.g. IVocks and ManrJ 120091 IBian and Kontarl 
UnSB lThrelfall, McClements, and de Moortel[ [5ji nilBian and Kontarl I50TT1 ITan et aU | 
I201ip . The simulations bv lReid and Kontarl (|2010p suggest that the suppression 
of Langmuir waves changes the electron beam transport. However, the detailed 
evolution of the electron spectra from the Sun to the Earth has not been per- 
formed. The question of electron transport and associated Langmuir waves in 
the inner heliosphere becomes particularly important in the view of anticipated 
observations by ESA's Solar Orbiter and NASA's Solar Probe Plus. 

In this paper, using a time-dependent injection, we investigate the detailed 
evolution of energetic electrons in the inner heliosphere. Specifically, we investi- 
gate the evolution of the spectral index change as the function of distance for 
both the electron fluence and the electron peak flux, the energy loss rate as a 
function of distance and the role of electron injection time during a solar flare. 



2. Electron transport of deka-keV electrons: model 

To model the transport of an electron beam from the Sun to the Earth we have 
taken a similar approach to the work that was done in lReid and Kontarl (|2010p . 
We use a set of two kinetic equations describing the motion of electrons along 
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the radially expanding magnetic field lines and their interaction with plasma via 
emission and absorption of Langmuir waves. Whilst the production of Langmuir 
waves is dealt with in a self-consistent manner, the non-linear effects of Langmuir 
waves scattering off ions and Langmuir wave coupling with ion sound waves 
is not considered (e.g. IKontar and Pecselii 12002 Zic bell et aUI20fTj) . Both the 
electrons and Langmuir waves exchange energy through the resonant interaction 
w pe = kv described by quasilinear terms [first term on the right hand sides 



(|Drummond and Pines, 1962} |Vedenov, Velikhov, and Sagdeev, 1962[)] 
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where f(v,r,t) is the electron beam distribution function, and W(v,r,t) the 
spectral energy density of Langmuir waves. The second term on the left had side 
of Equation ([I]) describes the propagation of electrons along the radial direction 
r. It also models a radial decrease in density of the electron beam as the beam 
fills the radial expanding magnetic field. This expansion is determined by the 
constant ro = 3.4 x 10 9 cm which is chosen to model a cone of expansion with 
an angle of 33.6°. The two terms on the right hand side of Equation ([1]) model 
the quasilinear diffusion of electrons in the presence of Langmuir waves and the 
collisional interaction of electrons off ions in the background plasma respectively. 

The second term on the left hand side of Equation describes the prop- 
agation of Langmuir waves along the radial direction r, where duj^jdk is the 
group velocity of Langmuir waves. The third term on the left hand side of 
Equation @ describes the refraction of Langmuir waves at the density gradients 
in the background heliospheric plasma. The first term on the right hand side of 
Equation ([2]) describes the growth rate of Langmuir waves from an unstable elec- 
tron beam through resonant interaction. j c and correspond to the collisional 
absorption of Langmuir waves and Landau damping respectively, where 7 C = 
7m e e 4 ln A/(mgi;f. e ) and = 2y / 7ru;p e (r) (v/vxe) 3 exp (— v 2 /v^ e ). The last term 
on the right hand side of Equation ([2]) represents the spontaneous wave genera- 
tion ( |Zheleznyakov and Zaitse~ 1970 Takakura and Shibahashi, 1976; Hannah, Kontar, and Sirenko, 2009] ) .| 

2.1. Density model in the solar corona and the inner heliosphcrc 



To model non-thermal electron transport through the heliosphere, taking into 
account self-induced Langmuir waves, an accurate model of the background 
electron density is required. Previous work ( Reid and Kontar, 2010 ) investigated 
the dependency of Langmuir wave generation from the Sun to the Earth with 
background density fluctuations. Such fluctuations are able to suppress Langmuir 
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wave growth and as such the turbulent nature of the solar wind must be taken 
into account for accurate modelling of electron beam propagation. Following the 
same approach as IReid and Kontarl (|2010p . we initially model the large scale 
decrease in electron density in the solar wind no(r) using the equations for a 



stationary spherical symmetric solution (Parker, 1958) with normalisation factor 
bv IMann etaU (fT9"9"9"|) . 

r 2 no(r)v(r) = C = const (3) 
#-lnf^=41n(^ + 4^-3 (4) 



where v c = v(r c ) = (ksT e / jimp) 1 / 2 , r c = GM s /2v 2 , T e is the electron temper- 
ature, M s is the mass of the Sun, m p is the proton mass and ft is the mean 
molecular weight. The constant appearing above is fixed by satellite measure- 
ments near the Earth's orbit (at r = 1 AU, n = 6.59 cm -3 ) and equates to 
6.3 x 10 34 s _1 . The density no(r) is calculated by numerical integration of Equa- 
tions ^ and (j3]). The density profile n(r) is then modelled by adding density 
fluctuations whose power density spectrum takes the form of a power-law with 
spectral index —5/3, 



n(r) = n (r) 



N 

l + C^Af/ 2 sm(27rr/A„ 

71 = 1 



(5) 



where A„ is the wavelength of the perturbations with random phase <p n , (3 = 5/3. 
C is a constant which defines the r.m.s. level of density fluctuations \J (An(r) 2 ) / (n(r)) 2 | 
or An(r)/n(r) for short. The range of wavelengths modelled is 10 7 < A„ < 
10 10 cm. We found previously (Reid and Kontar, 20101 that a high level of den- 



sity fluctuations (10%), an observed value near the Earth ( |Celnikier, Muschietti, and Goldman, 1987| ,[ 
all the way from the Sun to the Earth contradicts to the high level of Langmuir 
waves required to explain the observed solar type III radio bursts. The relative 
level of density fluctuations An(r)/n(r) is assumed to be decreasing towards the 
Sun from 10% at 1 AU using 

An(r) _ ( no(lAU) y AnjlAU) 
n(r) V Mr) ) n(lAU) [ ' 

where ip determines the rate at which levels of density fluctuations rise from the 
Sun to the Earth. It was found through a comparison with observations that a 
reasonable value for ip = 0.25. 



3. Initial electron injection 

To model the injection of energetic electrons into the solar corona, we use an 
electron distribution function f(v, r, t), which varies separately in velocity, space 
and time 

f(v,r,t) = g (v)h (r)i (t). (7) 
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The injected electron distribution function varies with velocity as a single power- 
law, consistent with HXR observations. 



n b (2S - 1) fv mhl \ 2S 
9o(v) = — (8) 



where 5 is the spectral index of the electron beam (in energy space), and rib is 
time integrated beam density taken as 10 7 cm~ 3 , such that realistic peak flux 
and fluence values are measured near the Earth ( |Lin, 1985] |Krucker et ai, 2 007; 



Krucker, Oakley, and Lin, 20091. w m i n is the minimum beam velocity which is 
set at 3 VTe, where vr e = 7.8 x 10 s cm s -1 is the background thermal electron 
velocity associated with a 2 MK plasma. u m in is defined close enough to vt £ such 
that Langmuir waves are suppressed via Landau damping. The highest velocity 
modelled is 2 x 10 10 cm s" 1 , which is equivalent to 115 keV under the non- 
rclativistic assumption. High velocities would require consideration of relativistic 
dynamics and are not required considering the range of energies usually involved 
in impulsive solar electron beam wave-particle interaction. 

The energetic electrons are injected in a finite volume near r = 0, which 
corresponds to height h = 5 x 10 9 cm (50 Mm) at a frequency of 415 MHz. The 
electron distribution takes the form of a Gaussian 



h (r) =cxp(^— J , (9) 

where d is the characteristic length scale of the electron acceleration region or 
one dimensional electron injection region, d is set to 10 9 cm (10 Mm) which is 
the longitudinal scale of the electron acceleration region found in the analysis 
by |Reid, Vilmer, a nd Kontar (|201ip . The length of the simulation box is set to 



1.2 AU which is a typical distance for electrons to reach the Earth, being larger 
than 1 AU due to the curvature of the assumed Parker spiral. 

Temporally, the electron distribution function also varies as a Gaussian with 
two characteristic timescales 

io (t) = :r irr7 j — ; cxp ( ~ {t ~ 2 t0)2 ) , (10) 

0.5 v / 7r(Tl + T 2 ) \ T / 

where r — t\ in the rise, for time t < to and r = t 2 in the decay time t > to- T\ 
and r 2 represent the characteristic rise and decay times with t\ < t 2 . By assum- 
ing a common acceleration region for both upward and downward propagating 
electron beams, observational values for n and t 2 can be obtained from hard X- 
ray measurements of solar flares ( |Holman et ai, 201 1| |Kontar et ai, 20TTj ) . We 
use Ti = 10 s and t 2 = 45 s which are the typical values for large flare hard X-ray 
rise and decay times. The value of to is set at to = 4ti to allow sufficient time for 
the rise of electron injection. The approximation of a Gaussian rise time will not 
capture the fine structure of electron acceleration in the corona, which is often 
evident in both X-ray and type III radio observations. However, it adequately 
reproduces the overall temporal behaviour of the electron acceleration. 

Observationally a correlation in the spectral index has been found for prompt 
electron events at 1 AU and HXR emitting electrons with energies > 50 keV 
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( |Krucker et al, 2007[ ). Such energies are likely to travel to 1 AU without al- 
tering the spectrum by generating Langmuir waves ( |Kontar and Reid, 20091 
|Reid and Kontar, 2010[ ). Using the same approximation of a common acceler- 
ation region, we set the injected spectral index S to be similar values found from 
hard X-ray observations. Although we note that a typical X-ray spectral indices 
show a so-called 'soft-hard-soft' behaviour in time (e.g. IParks and W incklcr, 
IT9"6T)1 iBenzl [T9T71 IHolman et oH l2"0"TTj) . where 6 starts high (soft spectrum), 
then goes low during the most intense period of emission (hard spectrum) and 
the returns high afterwards (soft spectrum). We vary S between 8 and 4 with 
the same time dependence used for electron injection (Equation (jlOp V 

We assume an initial spectral energy density of Langmuir waves at the thermal 
level 

viri + n\ k B T e uj pe (r) 2 ( v \ 

The thermal level of Langmuir waves represents the spontaneous emission of 
Langmuir waves from the background Maxwellian plasma. The spontaneous 
generation of Langmuir waves (the last term on the right hand side of Equation 
([2|)) only deals with the generation of Langmuir waves from the flare accelerated 
electron beam. Therefore to model the spontaneous emission of Langmuir waves 
at times t > we set the minimum level of Langmuir waves to be this thermal 
equilibrium W(v,r,t = 0). 



4. Energetic particle energy loss in the corona and in the heliosphere 

It has been found observationally that the peak flux and fluence of impul- 
sive electron beams can take the form of a broken power-law at the Earth 
(W ang, Fisk, and Lin, 1971] |Lin, 1974] ) . The break takes the form of a knee 
with lower energies having a smaller spectral index. A recent statistical study 
( |Krucker, Oakley, and Lin, 2009[ ) using the 3-D Plasma and Energetic Particle 



instrument (Lin et al., 19951 on the WIND spacecraft was carried out on 62 
impulsive events. The average break energy was found to be around 60 kcV 
although breaks were detected as low as 30 keV and sometimes above 100 kcV. 
Weaker events are also detected only at low and high energies depending upon 
the background flux of electrons. 

Recent work ( |Kontar and Reid, 2009| |Reid and Kontar, 2010[ ) explains that 
a broken power-law spectrum at the Earth is created from a single power-law 
injected electron spectrum if the generation and absorption of Langmuir waves 
are taken into account. Specifically the broken power-law forms as a combined 
process between the electron beam inducing Langmuir waves and the density 
inhomogeneity changing the phase velocity of the Langmuir waves. Above the 
break energy the electron beam is too dilute to generate Langmuir waves in 
the background plasma. However, below the break, Langmuir wave interaction 
with the electrons flattens the electron distribution function in velocity space. 
As the result, the peak flux and fluence spectrum flattens with the density 
inhomogeneity controlling the flattening. 
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An example of this behaviour is shown in Figure [T] at a distance of 21.6 R@ 
for a simulation using the equations from Section [5] and initial parameters from 
Section^ The top panel shows the lightcurves of the simulations (electron flux as 
a function of time) using the same energy bins as the WIND 3DP instrument. A 
power-law artificial background has also been added. One can observe higher en- 
ergy particles arriving before low energy particles. The middle and bottom panels 
show the electron flux and normalised spectral energy density as a function of 
electron energy at the different times. When the high energy particles (> 40 
kcV) arrive at 21.6 Rq we can observe no significant Langmuir wave growth. 
Consequently the high energy electron distribution is not changed through quasi- 
linear interaction. Lower energy particles (< 40 keV) have a higher flux and 
therefore are susceptible to wave-particle interactions. Note the broadened flux 
distribution in energy space at time t = 340 s and the corresponding high level of 
Langmuir waves above the thermal background. The peak flux spectrum (black 
line in the middle graph) has noticeably flattened below 40 keV. 

4.1. Collisional energy losses 

The collisional energy losses for the electron beam are modelled by the last 
term in Equation (JXJ) . Electron-electron collisions are important for low energy 
electrons in a dense background plasma of the corona near the injection site. As 
such, the electron beam only suffers significant energy loss from collisions near 
the acceleration site. Exactly how close to the acceleration site collisions play an 
important role is required for understanding the electron beam energetics. 

Figure [2] shows the ratio of electron beam energy with and without collisions 
for different energy bands as a function of distance from the acceleration site. 
We can see that collisions are really only important at small distances r < 1 
from the acceleration site. Moreover, the majority of energy losses occurs with 
low energy electrons < 6 keV. Electrons with energy E > 40 keV are largely 
unaffected by collisions. The resultant electron beam leaves the corona with 
approximately 0.7% of its injected energy. This percentage is heavily dependent 
upon the initial conditions of the electron beam, for example the starting height 
of the electron beam or the spectral index of the electron beam. 

4.2. Generation of Langmuir waves in uniform plasma 

To understand the energy loss an electron beam undergoes travelling through the 
turbulent medium of the solar wind we must first consider a uniform plasma. It 
was initially postulated ( jSturrock, 1964[ ) that an electron beam would lose all its 
energy to Langmuir waves over a very short distance. Later it was found through 
numerical simulations ( |Takakura and Shibahashi, 1976 ; Magel ssen and Smith, 1977] )| 
that an electron beam was able to propagate long distances by forming a beam- 
plasma structure. The structure travels through the heliosphere with the charac- 
teristic velocity of the energetic electrons. Langmuir waves, whose group velocity 
is orders of magnitude lower than the electron's kinetic velocity, are generated 
at the front of the electron beam and absorbed at the back of the electron beam. 
The generation of a beam-plasma structure is a result of the finite spatial width 
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Figure 1. The electron distribution function and Langmuir wave spectral energy density at 
21.5 solar radii from the injection point. Top: the electron beam flux as a function of time with 
energy channels denoted on the right hand side. The three black lines represent the times 140 
s, 210 s, 340 s used for the other two graphs. An artificial power-law background flux has been 
added for clarity. Middle: the electron beam flux as a function of energy at three points in 
time. The black curve represents the peak flux spectrum. Bottom: the spectral energy density 
of Langmuir waves normalised by the initial thermal level Wth a t the same points in time as 
the electron flux. Colours and symbols in the middle and bottom graph represent the same 
times. Note the lower energy electrons arriving at 21.5 solar radii at later times are associated 
with larger wave amplitudes. 



of an electron beam. The converse, a spatially uniform beam, would not undergo 
the bump-in-tail instability assuming a power-law electron distribution and the 
absence of collisions, as fast electrons would no longer be able to outpace slower 
electrons. 

It was mentioned as early as Takakura and Shibahashi ( 1976 ) that there is 
very little energy loss when an electron beam propagates in a beam-plasma 
structure. Indeed, the quasilinear terms of Equations ([T]) and @ conserve energy, 
momentum and density. Energy is temporarily exchanged between the particles 
and the waves. In the presence of a uniform background plasma all the energy 
that is transferred to Langmuir waves at the front of the beam should be re- 
absorbed by the electrons at the back of the beam. By running a simulation which 
just considered the quasilinear interaction of the electron beam we compared the 
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Figure 2. The electron beam energy as a function of distance from the acceleration site 
when electron-electron collisions are taken into account. The beam energy is normalised 
by the electron beam when no collisions are present. The red dashed line corresponds to 
the energy range -E m i n < E < 6 keV. The green dot dashed line corresponds to the en- 
ergy range 6 keV < E < 40 keV. The blue solid line corresponds to the energy range 
40 keV < E < 100 keV. 



energetics to a free-streaming simulation. There was no significant net loss of 
energy when waves are produced (under the assumption of no binary Coulomb 
collisions) for energies above where Landau damping is significant. However, 
the inner heliosphere is highly anisotropic so this simple situation where beam 
energy is conserved does not apply to type III producing electron beams. 

4.3. Generation of Langmuir waves and collisional damping of waves 

One process which can reduce the energy of a beam-plasma structure is collisional 
damping of Langmuir waves in the background plasma where the damping rate 

4 

is 7c ~ "2^) In A. Similar to particle collisions, the collisional damping of waves 
is heavily dependent upon the background density of electrons. As a result it 
has the greatest effect close to the Sun. However, an important difference is that 
collisional damping of waves is independent on the phase velocity of Langmuir 
waves. 

The energy contained in Langmuir waves is much smaller than the kinetic 
energy of the electrons. Therefore, we do not expect wave collisions to reduce the 
total energy in the beam-plasma structure to a similar extent as particle-particle 
collisions. Given the initial conditions described in Section [3] we find that the 
effect of wave collisions is negligible on the total energy. Wave-particle instability 
in the beam of particles does not occur instantly but takes around 2 Rq of 
propagation to be noticeable (see Section [7] for a discussion). The background 
density at such distances is too small for significant collisional absorption of 
Langmuir wave energy. The collisional damping of waves will play an important 
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role deep in the corona. However, it is not the dominant energy loss process for 
electrons travelling through the heliosphere. 

4.4. Energy losses due to inhomogeneous plasma 

The spatial gradient of the background electron density plays an important role 
in the dynamics of Langmuir waves. It has been shown before (e.g. IKontarl 
200la|) that Langmuir waves arc refracted to different k-vectors by the inho- 
mogeneous background electron density gradient (Term 2 of the left hand side 
of Equation ©). The radial dependence in this process is mostly governed by 
the characteristic scale of plasma inhomogencity L(r) = 2n e (r){dn e {r) / dr)^ 1 . 
For a simple solar wind plasma model where density only decreases L(r) is 
strictly negative. When fluctuations are added to mimic the turbulent nature 
of the solar wind L(r) not only becomes positive in parts but varies more in 
magnitude. Such fluctuations in the background electron density can suppress 
Langmuir wave growth by moving Langmuir waves to higher and lower k-vectors, 
out of resonance with inducing electrons. 

To examine how the electron beam varies with distance the total beam-plasma 
energy has been plotted for various energy ranges in Figure [3J Values were nor- 
malised using the simulations from Section |4~21 where wave interaction was taken 
into account but not the role of inhomogeneous plasma. We define three energy 
channels to analyse the energetics E m - m < E < 6 keV, 6 keV < E < 40 keV 
and 6 keV < E < 40 keV. Resonant Langmuir waves whose phase velocity is the 
same as the electron velocities within the three energy ranges are included in 
the energetics. 

In Figure [3] we observe no net change in energy below 6 keV. The bulk of en- 
ergy is stored in electrons near E m - m where Langmuir waves are absorbed due to 
Landau damping by the background Maxwellian plasma. Between 6 and 40 keV, 
beam-plasma energy starts to decrease around 2 Rq when significant levels of 
Langmuir waves become induced by the electron beam. By 20 Rq this energy 
range has lost approximately half the energy it would have without the effect of 
inhomogeneous plasma. Conversely the higher energy channel between 40 and 
100 keV gains energy as Langmuir waves are absorbed. This occurs at a greater 
distance than 2 Rq because it is linked to positive values of L(r) which moves 
wave energy to lower values in k-space (higher phase velocity). It is not until the 
beam propagates further away from the Sun that effects from the background 
electron density turbulence are significant over the radially decreasing density 
decrease (see lReid and Kontar! I201dl for further details). When positive values of 
L(r) become significant, wave energy is moved to higher phase velocities and can 
be absorbed by higher energy electrons. Positive values of L(r) can be observed 
as the spikes in the high energy channel and correspond to waves at high phase 
velocities. Thus, inhomogencity causes a spectral shift in energy. 

4.5. Energy loss from the Sun to the Earth 

After looking at the energy loss from the individual terms, we now consider prop- 
agation from the Sun to the Earth using all physical processes defined previously 
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Figure 3. The total energy in the beam-plasma structure (waves and particles) as a function 
of distance when background plasma inhomogencity is taken into account. Energy is normalised 
using the simulation where Langmuir waves are induced from the quasilincar term but the inho- 
mogeneity term was ignored. The red line corresponds to the energy range i? m i n < E < 6 keV. 
The green line corresponds to the energy range 6 keV < E < 40 keV. The blue line corresponds 
to the energy range 40 keV < E < 100 keV. 

in Equations ([I} and ©. Figure 0] shows how the energetics in the previously 
defined energy channels varies with radial distance from the Sun normalised 
by the case where only free-streaming electrons are considered. We observe the 
same initial high level of energy loss of the electron beam in the corona which 
was predicted in our analysis of the coulomb collisions. The dense plasma of the 
corona damps the low velocity electrons which contain the majority of the total 
energy in the electron beam. 

After 1 Rq the background plasma rarefies and inhomogeneity plays the 
dominant role for the electron beam energy loss. By 1 AU the energy range 
6 keV < E < 40 keV has lost nearly 1 order of magnitude of energy due to 
inhomogeneity with respect to the free-streaming case. Conversely the higher 
energy range 40 keV < E < 100 keV has gained a small amount of energy by 
absorbing wave energy refracted to higher phase velocity. 

5. Electron spectrum evolution 

Although a break in the power-law spectra appears to be formed in the electron 
spectrum from propagation effects, it is not obvious how the break forms as a 
function of distance from the Sun. Our previous work ( |Reid and Kontar, 2010| ) 
suggests that high energy (> 20 keV) electrons only participate in wave- particle 
interactions close to the Sun before their flux rarefies due to the expansion 
of the magnetic field. In this scenario the power-law break would be formed 
quite early in the propagation of the electron beam towards the Earth. Any 
predictions concerning electron beam dynamics between the Sun and the Earth 




SOLA: inner_helio_decl7.tex; 25 September 2012; 0:33; p. 11 



Rcid and Kontar 




0.1 1.0 10.0 100.0 

Distance from Injection Site [R ] 



Figure 4. Total energy of the beam-plasma structure (waves and particles) as a function 
of distance from the Sun in three different energy bands. The energy is normalised using 
the simulation where electrons are free-streaming only. The red line corresponds to the energy 
range E m - ln < E < 6 keV. The green line corresponds to the energy range 6 keV < E < 40 keV. 
The blue line corresponds to the energy range 40 keV < E < 100 keV. 

is particularly relevant with the upcoming missions of Solar Orbiter and Solar 
Probe Plus later this decade. As such we investigate how the electron spectrum 
both in peak flux and fluence (time integrated flux) varies radially from the Sun 
to the Earth in our simulations (Figure [5]). 

To investigate the evolution of the electron beam spectrum, we fit both the 
peak flux and fluence with a broken power-law as a function of distance. The free 
parameters for the fit are the break energy, the flux or fluence at the break energy 
and the spectral index both below and above the break energy. The fit is made 
between 6 and 100 keV. Figure [5] shows the variation of three free parameters 
for the fit as a function of distance. The weighting for the fit was 1/y 2 such 
that x 2 = J2(v ~ Vfit) 2 /y 2 allowing the closest logarithmic fit to the data. By 
assuming an implicitly good fit we have scaled the one sigma errors with the 
chi-squared values using a c = ay/ x 2 /DOF, where DOF represents the degrees 
of freedom in the fit, discerned by the number of points minus the number of fit 
parameters. To analyse Figure |5] we will consider both the peak flux and fluence 
at specific distances from the Sun that arc plotted in Figure [5] 

Initially, as shown in the previous section on energetics, the collisions domi- 
nate the transport of the low energy electrons inside 1 Rq. Such a scenario will 
not be ubiquitous of all electron beams. However, due to the high characteristic 
time of the temporal electron beam injection and the starting background elec- 
tron density we chose it is expected. Figure [5] shows a snapshot of the peak flux 
and fluence spectrum at 1.5 Rq. The data has been logarithmically binned in 
energy space to show a clear representation. The green dashed line shows a best 
fit to the data using a broken power-law function. The break energy is very low 
(< 20 keV) on account of the spectral change from Coulomb collisions affecting 
low energy electrons. The spectral index above the break energy for both the 
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Figure 5. Fluence (red triangles) and peak flux (blue diamonds) spectrum of energetic elec- 
trons at different distances from the injection site. The bottom right panel shows the spectrum 
at 1.2 AU. A green best fit line has been drawn using a broken power-law function for all 
curves. 



peak flux and fluence is around 4.5 which is exactly what would be expected for 
free streaming electrons with initial injection distribution described by Equation 

We have already shown that wave-particle interactions have become ener- 
getically important for the electron beam by 2.5 Rq. Again, we highlight this 
distance is defined by the initial injection parameters of the electron beam sim- 
ulated. When wave-particle interactions become important the spectral index 
of the beam starts to change below 40 keV. Figure [5] shows a snapshot of the 
spectrum at 7.5 Rq. A bump in the fluence spectrum can be seen around 40 keV 
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Figure 6. The broken power-law fit parameters for the fluence (red) and peak flux (blue) 
spectra as a function of distance from the Sun. Top: spectral index below the break energy. 
Middle: spectral index above the break energy. Bottom: break energy. 

along with spectral flattening below 30 keV. The peak flux also shows spectral 
flattening below 40 keV. The spectral flattening is enough that the reduction 
in flux from collisions at low energies < 10 keV is no longer visible. The bump 
in fluence spectrum is caused by the absorption of wave energy which has been 
shifted to higher phase velocities by positive background electron density gradi- 
ents. The bump has caused an increase in the fluence spectral index above the 
break energy from 4.5 to around 5. The break energy for both the peak flux and 
fluence has increased to between 30 and 40 keV. 

From Figure [5] we observe the spectral index below the break energy de- 
creasing for both peak flux and fluence spectra till about 20 Rq. We note that 
the energy spectrum observed in the simulations is not an exact broken power- 
law but more of a triple power-law. However, it can still be approximated with a 
double power-law fit. This can be seen from the poorness of the fit represented by 
the increased errors presented in Figure |5] After 20 Rq the spectral parameters 
change slowly. The peak flux and fluence spectra in Figure [5] are very similar at 
distances of 50 R Q (0.23 AU), 150 R Q (0.7 AU) and 258.6 R (1.2 AU). The 
spectral index below the break for the peak flux decreases the most from around 2 
to 1.2. while the fluence spectral index below the break remains almost constant. 
The density inhomogeneity also leads to re-acceleration of beam electrons evident 
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Figure 7. The peak flux and fluence of the electron beam at 1.2 AU with density fluctuations 
(filled symbols) and without fluctuations (empty symbols) in the heliosphere assumed. Without 
considering density fluctuations the peak flux has a flat distribution between 10 and 40 keV 
that does not agree with observations near the Earth. 



in the bump in the fluence spectrum. This is related to Langmuir waves shifted 
to smaller k (larger velocities) by the density inhomogeneities and consequently 
reabsorbed by the higher energy electrons, which results in appearance of acceler- 
ated electrons ( jNishikawa and Ryutov, 1976| |Kontar, 2001b[ ) . The bump in the 
fluence also makes the break energy in the fluence spectrum slightly ambiguous 
compared to the sharp break in the peak flux spectrum. 



6. Electron density turbulence 

The level of density fluctuations An(r)/n(r) directly affects the amount of 
Langmuir waves produced by the electron beam as it travels through the inner 
heliosphere. If no density fluctuations are present a higher level of Langmuir 
waves is produced. To highlight their importance we ran a simulation which only 
modelled a radially decreasing density and compared the spectrum at the Earth 
to a simulation which included density fluctuations. The corresponding peak 
flux and fluence spectrum is presented in Figure [7J where the filled and empty 
symbols represent the case with and without density fluctuations respectively. 
The increased level of Langmuir waves from a lack of density fluctuations caused 
increased flattening of the spectrum below the break energy in both fluence and 
peak flux. Indeed the flattening is so extreme that the peak flux spectrum is 
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a plateau (spectral index around 0) between the energies 10-40 keV. Such a 
situation is never observed in-situ with spacecraft. 

The cause of such a flat spectrum is the presence of an anomalously large 
magnitude of Langmuir waves when the electrons below 40 keV arrive at the 
Earth. A combination between a plateau forming in the electron distribution 
function and wave refraction causes this constant maximum between 10 and 40 
keV. The observed increase at 10 keV is caused by Landau damping absorbing 
Langmuir waves at low energies. The electron distribution function consequently 
cannot spread as much in energy space and, as number density is conserved, we 
observe both as increase in peak flux and flucncc. Whilst this scenario of zero 
density fluctuations is unrealistic for the solar wind, a relatively low level of 
density fluctuations could be responsible for observations of solar beam-plasma 
structures associated with type III radio bursts at distances as large as 4.3 AU 
(|Buttighoffer et ai, 19951. 



7. Discussion and conclusions 

We have investigated the changing energetics of a solar electron beam as it travels 
through the inner hcliosphere. The majority of energy loss for an electron beam is 
shown to happen in the dense solar corona and comes from Coulomb collisions 
of electrons with ions. When the electron beam propagates into the rarefied 
inner heliosphere, provided there are enough electrons, wave-particle interaction 
becomes the energetically dominant effect. Specifically for our simulations the 
energy range 6 keV < E < 40 keV lost a large percentage of energy due to 
resonant interaction with Langmuir waves and the subsequent Langmuir wave in- 
teraction with the inhomogeneous background plasma. The loss of energy caused 
a spectral flattening of the initial power-law electron distribution which resulted 
in a broken power-law (with break around 40 keV) being a very good fit to the 
data. We have shown how the broken power-law fit parameters vary as a function 
of distance, finding the best fit to the data around 1 AU and the worst fit around 
10-30 Rq. Finally, we have shown that modelling the turbulent fluctuations in 
the inner heliosphere is essential to obtain both an electron distribution and a 
level of Langmuir waves that is consistent with observations near the Earth. 

It should be stressed that wave-particle interaction not occurring above 40 keV 
is entirely dependent upon the initial beam characteristics. Increasing the beam 
density or decreasing the injected spectral index, characteristic time or size will 
all result in a higher number density of electrons that can become unstable 
to wave-particle interactions and can consequently have a higher break energy. 
Although the parameters we have chosen do reflect typical flare parameters 
(see Section [3]) we could have run the simulations in a slightly different area 
of parameter space which would have produced higher or lower break energies. 
Other ways to increase the break energy is to inject the electron beam at lower 
background electron densities or reduce the expansion factor of the magnetic 
field. 

We showed how the electron beam had to travel around 2 R Q before it 
would become unstable to the generation of Langmuir waves. Previous work 
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( |Reid, Vilmer, and Kontar, 201 1[ ) showed that the instability distance of an elec- 
tron beam was heavily dependent upon the injected spectral index and the 
longitudinal extent of the acceleration site. However, the modelled injection of 
the electron beam was instantaneous in time. In this paper, we have considered 
electron beam injection which is time dependent, with a Gaussian profile and 
characteristic times, described by Equation (|10j) . Using a similar analysis to 



Rcid , Vilmer, and Kontar ([2011[) . we can analyse how a temporal injection alters 
the instability distance of an electron beam. If we assume an injected distribution 
of the form 

f{v,r,t) = g (v)exp(-\r\/d)(TTT)-i exp(-\t - t \/r). (12) 

where both the spatial and temporal form of the electrons is an exponential 
distribution function. We obtain an instability criteria of the form 

htypein = {d + vr)2S + h acc (13) 

where h acc is the injection height from the photosphere and h ty peiii is the height 
that a large level of Langmuir waves are induced and type III emission is ob- 
served. The instability distance is now dependent upon the spectral index 5 (26 
is the spectral index in velocity space) and the highest of the terms in d + vt. 
Assuming an injection site of 10 9 cm, spectral index S = 4, a characteristic 
velocity of 5 x 10 9 cm s _1 and a injection time of 10 seconds we find that Lang- 
muir wave generation will not happen for around 5.5 Rq. The exact distribution 
function used in the simulations varies slightly from Equation (|12[). but we find 
Langmuir wave generation occurring around 2.5 Rq when we consider quasilincar 
interaction only. When we consider the collisional term, we find Langmuir wave 
growth develops much quicker. The collisions cause low velocity electrons near 
VTe to lose energy, creating a positive gradient in velocity space earlier than 
time-of- flight would otherwise produce ( |Hannah and Kontar, 20lTj ). 

Finally we note that the spectrum evolves quicker at the distances closer to 
the Sun, so that future inner heliosphere missions can diagnose the evolution of 
energetic electrons, better restricting the possible parameters of the plasma and 
energetic electrons. For the parameters adopted in the model, the simulations 
predict weakly changing spectral indices and break energy after around r ~ 
20-R©, and the decrease of the spectral index below the break at the distances 
r < 2Oi?0. The spectral break, well pronounced at 1 AU, is less evident at the 
distances r < 2Oi?0, so the peak flux spectrum will be closer to the electron flux 
spectrum inferred from hard X-ray observations. 
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